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We analyze the creation of spin squeezed atomic ensembles by simultaneous dispersive interactions 
with several optical frequencies. A judicious choice of optical parameters enables optimization 
of an interferometric detection scheme that suppresses inhomogeneous light shifts and keeps the 
interferometer operating in a balanced mode that minimizes technical noise. We show that when 
the atoms interact with two-frequency light tuned to cycling transitions the degree of spin squeezing 
scales as ~ 1/d where d is the resonant optical depth of the ensemble. In real alkali atoms 
there are loss channels and the scaling may be closer to ~ l/\/d. Nevertheless the use of two- 
frequencies provides a significant improvement in the degree of squeezing attainable as we show by 
quantitative analysis of non-resonant probing on the Cs Dl line. Two alternative configurations 
are analyzed: a Mach-Zehnder interferometer that uses spatial interference, and an interaction with 
multi-frequency amplitude modulated light that does not require a spatial interferometer. 

PACS numbers: 42.50.Nn, 42.50.Lc, 32.80.Qk, 03.65.Ta 



I. INTRODUCTION 

Coupling between light beams and atomic ensembles 
is of interest for processing and storing quantum infor- 
mation, and for enabling high precision measurements 
of fundamental physical quantities [ij. Recent develop- 
ments in atomic clocks have demonstrated a measure- 
ment uncertainty that is limited by the quantum projec- 
tion noise of atomic spin measurements |2|. Reductions 
in the measurement uncertainty may be achieved by us- 
ing spin squeezed states (SSS) of atomic ensembles d, 3 ■ 
Such states were generated using an off-resonant quan- 
tum nondemolition (QND) interaction with a coherent 
light beam[^, Q. Preparation of an atomic sample in a 
SSS reduces the variance of a projective measurement of 
the spin by a factor of = 1/(1 -|- /t^) below that of an 
ensemble prepared in a coherent spin state (CSS). Here 
k; is a constant proportional to the light matter interac- 
tion strength. Development of techniques for generating 
strongly squeezed atomic samples is therefore of great 
interest as a route to improving the precision of atomic 
clocks. 

Spin squeezed states can be generated via a QND inter- 
action described by a Hamiltonian of the form TYqnd ^ 
FzSz, with F and S referring to spin degrees of freedom 
of the atoms and the light respectively. In this paper 
we will focus on the situation where the atomic pseu- 
dospin F is defined in the basis |3) = |/ = 3,m/ = 0) 
and |4) = |/ = 4, = 0) corresponding to the clock 
transition between / = 3 and / = 4 ground state hy- 
perfine levels in Cs. The optical pseudospin operator S 
may be defined in a basis of polarization, spatial, or fre- 
quency modes of the light. One of the challenges en- 
countered in preparation of spin squeezing is the fact 
that the above QND Hamiltonian is only an approxima- 



tion that neglects additional aspects of the light-matter 
interaction which serve to reduce the usable amount of 
squeezing. For example in the case of polarization depen- 
dent optical Faraday rotation there are nonlinear terms 
in the atomic tensor polarizability which give a non QND 
like interaction^, [l] leading to decoherence of atomic su- 
perposition states. We will consider z polarized light 
beams as was used in our recent observation of Rabi 
oscillations on the Cs clock transition^. The z polar- 
ized light couples to the atomic basis states without any 
nonlinear terms. Nevertheless spatial inhomogeneity of 
the light-atom coupling strength leads to inhomogeneous 
atomic phase shifts. It is still possible to obtain a strongly 
squeezed ensemble characterized by a nonsymmetric en- 
tanglement measure in this situation [loj. However, the 
presence of inhomogeneous coupling is problematic in the 
context of reducing projection noise in atomic clock ex- 
periments, since it limits the fidelity with which ensemble 
rotation operations can be performed 11 1 . 



In this paper we study the use of multi-frequency 
light beams for creating SSS. It was first shown in |12l | 
that also with multiple probe frequencies a QND in- 
teraction can be obtained. Here we consider a pair of 
two-frequency amplitude modulated light fields that are 
analogous to a carrier suppressed frequency modulation 
(FM) spectroscopy [Tsj . We use carrier frequencies sym- 
metrically placed with respect to an atomic resonance 
to engineer an effective QND interaction, while canceling 
inhomogeneous light shifts. In the spirit of Ref. [ic| mul- 
tiple frequencies encoded in a light beam with a common 
spatial mode provides a convenient method to ensure the 
matched interactions needed for nonsymmetric entangle- 
ment generation. We proceed in Sec. |TT] by recalling 
the form of the interaction between an atomic ensemble 
and an off-resonant light field and estimate the degree of 
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FIG. 1: (color online) Interferometric setup for QND mea- 
surement of atomic spin on the Cs Dl line. The coherent 
state input mode c is split into two equal intensity parts at 
the first beamsplitter, and then recombined at the output 
beamsplitter to give output modes di, d2 that are measured 
by photodetectors to generate a difference current i_. 



squeezing obtainable by probing of the Cs Dl line. We 
show that the interaction with a linearly polarized single 
frequency probe beam suffers from inhomogeneous light 
shifts. The inhomogeneous shifts can be eliminated us- 
ing two probe beams of different frequencies in a Mach- 
Zchnder interferometer as discussed in Sec. IIIII Pro- 
vided the number of photons is large compared to the 
number of atoms this interaction has the potential for 
producing SSS without unwanted inhomogeneous light 
shifts. For an idealized light-atom interaction the two- 
frequency technique leads to spin squeezing that scales 
as ~ 1 /d where d is the resonant optical depth of the 
ensemble. Although loss channels on the Cs Dl line limit 
the asymptotic squeezing to ~ 1/Vd we show never- 
theless that the quantitative performance is better than 
for one-frequency probing. 

In Sec. IIVI we present an alternative configuration 
which eliminates the spatial Mach-Zehnder interferom- 
eter in favor of frequency domain phase shifts. This is 
attractive since it removes the requirement of mechani- 
cal stability inherent in using an interferometer. We con- 
clude in Sec. Ivlwith a discussion of the results obtained. 



II. QND MEASUREMENT WITH A SINGLE 
PROBE BEAM 



Consider an interferometric measurement of the collec- 
tive atomic spin using the arrangement shown in Fig. [T] 
This type of setup has been analyzed in detail in several 
papers (l3.[l5j and we will limit ourselves to a brief discus- 
sion in order to define notation to be used in what follows. 
An input beam in a coherent state is split into two parts 
with equal amplitudes. A cloud of atoms placed in one 
arm of the Mach-Zehnder imparts a phase shift on the 
light. The phase shift is transformed into an electronic 
signal by measuring the difference of the photocurrents of 
the output detectors. Essentially this configuration was 
used in Ref. Q except that the probing light was tuned 
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FIG. 2: (color online) Refractive index and differential light 
shift near the Dl line for z polarized light, an atomic density 
of 1 X lO" cm~^ prepared with (Fz) = 0, A = nw'^ with 
It) = 20 /im, and an optical power of P = 10"® W. 



close to the Cs D2 line, instead of the Dl line considered 
here. 

The QND interaction generating spin squeezing in 
our system is governed by the Hamiltonian TYqnd ^ 

F4Nph/2 + S,) where = Eti and & are 
operators for the z component of collective atomic fl6j 
and photonic spins respectively and Np}^ is the pho- 
ton number operator (for convenience we set h = 1). 
The dimensionless single particle pseudospin operators 



are fi = 
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2 \ '"44 '"ss / ^'^^ ^'^'^ atoms, while the contin- 
uous operators describing the light field are written as 
5^ = i 4- fo^a^ tp, Sy = {^'^b ~ t^, = 
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ments of the single atom slowly varying density operator, 
d is the annihilation operator for the field that interacts 
with the atoms, and b is the annihilation operator for the 
local oscillator field in the lower arm of the interferome- 
ter which has no direct interaction with the atoms. The 
above definitions are supplemented by the number oper- 
ators for atoms and photons: A^at = 



P\i 



and TVph = <p ( d^d + b^b 
light pulse. 



with tp the duration of the 
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The probe beam refractive index and differential light 
shift of the clock states |3), |4) are shown in Fig. [5] for 
light of wavelength A and frequency lo that is near res- 
onant with the Dl line of Cs. A probe beam tuned to 
the zero index shift point indicated in Fig. [2] receives 
equal and opposite phase shifts from population in the 
/ = 3 and / = 4 states. The effective refractive index 
is therefore proportional to the population difference of 
these states and is given by = 1 — ^2{Fz)k. When 
the probe light is tt polarized (linearly polarized along z) 
the interaction constant is 
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A is the transverse area of the light beam and atomic 
ensemble, I is the length of the ensemble, Ajk = ui — 
(■^6pi/2,/'=j ~ -£^661/2 J=fe)/^ is the detuning of the probe 
beam from the corresponding optical transition, and 7 
is the radiative linewidth (FWHM) of the excited states. 
We assume that the atoms are cold so that it is not nec- 
essary to account for the presence of Doppler broadening. 

In the limit of large detuning where the photon scat- 
tering rate and accompanying population changes are 
small the pseudospin operators transform as F°"' = 
R,(^at)F•^ S°-' = R,(^ph)S•^ with 
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It follows from the form of TYqnd that the rotation angles 
are 
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FIG. 3: (color online) Bloch sphere representation of CSS 
of atoms and light. The top two rows show the initial states 
and their fluctuations. The output states after the interaction 
show that the atoms and light are rotated about the 2 axis 
leading to unequal fluctuations transverse to the mean spin 
directions. Measurement of the light with quantum limited 
uncertainty projects the atoms into a SSS with the variance 
of Fz reduced by a factor of 



0,t = 2K(iVph/2 + 5;") (2a) 
^ph = -2kFl'\ (2b) 

These angles characterize the strength of the light-atom 
coupling. 

Consider atoms and photons prepared in CSS's as 
shown in Fig. [31 The atomic pseudospin is aligned such 
that (Fi") = (Fi") = 0, and (F'") = {N,t)/2 = N,t/2. 
For the atoms, we may assume that we prepared the CSS 
by starting with a definite number of atoms in |3) so that 
initially {Fz) = —Nat/2 and then used a perfect 7r/2 pulse 
to create the state with (Fy) = Nat/2. The variances of 
the prepared state are {{AF'^'f) = ((AFj")^) = iVat/4 
and ((AF^")2) = 0. As for the Hght, the input port of 
the interferometer divides the light equally between the 
two arms giving (^j,") = (^f) = 0, 01^) = (7Vph)/2 = 
7Vph/2 and input variances ((AS'i")^) = ((ASj,")^) = 
((AS*"')^) ~ A^ph/4. With these initial conditions we have 
(^ph) = but (^at) 7^ due to the presence of a nonzero 
differential light shift of the atomic states. We will return 
to the significance of the light shift below. 

To lowest order in the interaction strength we find the 



output variance of the light is 

{{AS;^'f) - {{AS'j^f) + {2kf{A{FrSTf) 

= ^{l + ^). (3) 

where = jk^N^tNph- 

To be of use, the operators for the light interacting 
with the atoms in one arm of a Mach - Zehnder inter- 
ferometer as shown in Fig. [H must be related to a light 
observable at the interferometer output. To that end, we 
choose the phase of the interferometer so that the pow- 
ers in the two outputs are equal. At this position the 
interferometer has the greatest sensitivity to anticipated 
phase changes from the atomic sample and, additionally, 
any classical noise on the laser power is rejected. The 
field operators at the output beam splitter shown in Fig. 
[T] are related to the internal fields by di — "^(o- + 
^2 = "71^^ ~ ^Consequently, the components of the 
light operator S transform as Sdx = Sz, Sdy = —S^, 
Sdz = ^Sy. where subscript d refers to the interferom- 
eter output fields. The output beams are detected and 
the photocurrents subtracted. The difference photocur- 
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rent is proportional to L 
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Combining with Eq. ([3]) we have = and 
((Az_)2) = 4((A^™*)2) = iVph(l + 



2Sy. 



(4) 



When there are no atoms the variance of the measured 
difference current is given by iVph, the coherent state 
rcsuh as expected. When atoms are present (k^ > 0) 
the variance increases Unearly with the number of atoms, 
which is just the projection noise of a CSS. 

A single quantum limited measurement of the differ- 
ence photocurrent with variance Npi^ represents a reduc- 
tion by a factor of 1 -I- compared to the variance given 
by Eq. As is shown pictorially in Fig. [3] the mea- 

surement projects the atoms into a spin squeezed state 
(SSS) with the variance of the z component reduced by 
the same factor, such that 
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The SSS is characterized by the degree of squeezing [3, 0] 
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A detailed discussion of the projective reduction of the 
atomic variance using a wavcfunction formalism can be 
found in In order to reduce the uncertainty in a 

Ramsey measurement of an atomic clock frequency addi- 
tional operations are needed which include interchanging 
the variances of Fx and Fz as described in Ref. [isj . 

The degree of spin squeezing given by Eq. ([6]) neglects 
the deleterious effects of inelastic light scattering which 
reduces the magnitude of the coherent spin state and 
adds noise to the z components of the pseudo spins. It is 
well known [17| that the maximum attainable spin squeez- 
ing accounting for light scattering scales as ~ 1/Vd 
with d the resonant optical depth of the atomic sam- 
ple. An exact calculation of the degree of spin squeezing 
in a real atomic system is very cumbersome. Previous 
work has provided analytical results in a Gaussian ap- 
proximation supplemented by numerical analysis to ac- 
count for atomic decay and light scattering [ij, [11]. The 
Gaussian state based calculations assume an idealized 
two-level atomic structure. The effect of atomic redistri- 
bution to other internal states due to light scattering was 
accounted for in [l^ for the case of ^''Rb probed on the 
D2 line using an approximate analysis valid for not too 
large decay rates. Here we follow the spirit of [Tol in the 
setting of the Cs Dl line. 

With TT polarized light tuned between the resonances 
as shown in Fig. [5] it is readily shown that ~ dri/2 
where d = aoN^^t/A, ctq = A^/27r, and rj is the integrated 
probability that an atom scatters a photon during the 
probing pulse. The result of the coherent QND interac- 
tion can thus be written as .^^ = 1/(1 -|- ■^drj). Inelastic 
scattering events result in either decoherencc and return 
of an atom to its original state with probability r^dc , or 
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FIG. 4: Optimized spin squeezing and scattering parameter 
as a function of optical depth d from Eq. (0 ■ The inset shows 
the variation of ^'^ with rj at d = 100. Realistic experimental 
parameters for Cs atoms in an optical dipole trap as in[^ are 
f5) A = 0.894 fim, A = vrr^ r ^ 50 fim, I = 2{Tvr^/\) = 1.8 cm, 
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22. With iVph = 9.5 x 10'" we get ?? = 0.15 



loss to states with m/ = ±1 with probability rji. These 
"loss" states couple to the probing light with slightly dif- 
ferent strengths. For the Cs Dl line and tt polarized light 
the coupling is about 6% weaker for my = ±1 than for 
mf ~ 0. As we have defined the atomic pseudospin in 
the basis oi nif =0 states we will consider population 
of mf 7^ states as a loss mechanism and ignore the 
coupling of these states to the probing light. In practice 
this assumption may be made realistic by interspersing 
the QND interaction with cleaning steps that remove the 
population of nif ^ states. This could be done by, 
for example, coherently shelving the populations of the 
basis states to other levels and blowing away any popu- 
lation in the m 7^ states using resonant light. Account- 
ing for the relevant Clebsch-Gordan coefficients we find 
r// = 277/3, rydc =^7/3 and a short calculation then gives 
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For 3> 1 the squeezing is optimized for 770 ~ y^3/2d 
which gives ^^j^^ ~ y^32/3d. Figure [4] shows the degree 
of spin squeezing as a function of optical depth and scat- 
tering probability. We see that at d = 100 the optimum 
is to set 7^0 — 0.10 which gives ^^j^^ ~ 0.30. 

The above discussion is still highly idealized in that 
it assumes a uniform interaction strength for all atoms 
in the sample. The zero phase shift frequency shown in 
Fig. [2] imparts unequal light shifts to the clock states. 
In a practical situation with a probing beam of Gaussian 
profile the strength of the light field will vary across the 
sample leading to inhomogeneous broadening and rapid 
loss of coherence between the clock states. The resulting 
dcphasing of Rabi oscillations on the clock transition has 
been shown to be well described by a model that accounts 



5 



for a Gaussian beam profile and a Gaussian distribution 
of atoms in the ensemble ll| . In order to eliminate the 
inhomogeneous broadening it is necessary to use a probe 
frequency that results in equal light shifts for both clock 
states. This is possible using tt polarization on the D2 
line or linear polarization at an angle of 45° from z on the 
Dl line [3 • Unfortunately, as shown in Q. the frequen- 
cies for which the light shifts are equalized arc relatively 
close to resonance leading to strong photon scattering, 
and correspond to a nonzero phase shift of the light, so 
the interferometer operates in an undesired unbalanced 
configuration. 



III. QND MEASUREMENT WITH TWO PROBE 
BEAMS 



We now show that it is possible to eliminate the in- 
homogeneous light shifts, and operate the interferome- 
ter in a balanced configuration, while retaining the free- 
dom of choosing the detuning to optimize the interaction 
strength. To achieve this we use two tt polarized beams, 
one of frequency tuned close to the / = 3 /' = 4 
transition and one of frequency 074 tuned close to the 
/ = 4 ^ /' = 3 transition. We introduce two sets of 
continuous operators S^x, S^y, S^z and Sax, S^y, Saz for 
the two light fields respectively. These are defined in the 
same way as in the previous section with the replace- 
ments a — > 0,3, 6 — > 63, etc. . 

The light atom interaction is now characterized by four 
coupling constants corresponding to the interaction of 
each frequency with each of the ground states. We will 
be interested in detunings such that the interaction of W3 
light with population in / = 4, and the interaction of 
uja light with population in / = 3 is about one hundred 
times weaker than the interaction of each frequency with 
the population of the near resonant levels. We therefore 
only need to consider the two coupling constants to the 
near resonant levels 
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where A, 



Two conflicting requirements arise when using multi- 
ple probe beams. On the one hand we would like the 
total differential light shift to vanish. As seen in Fig. [2] 
this implies that A34 and A43 should have the same sign. 
On the other hand, in order to main to maintain equal 
intensities at both output ports of the interferometer, 
which minimizes technical noise in the detection process, 
the two probe beams should acquire opposite phase shifts 
from their near resonant level implying that A34 and A43 
should have opposite signs. We can solve this difhculty 
by choosing A34 , A43 to have the same sign but let them 
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FIG. 5: (color online) Two frequency QND measurement of 
atomic spin. The probe frequencies are injected into different 
input ports as discussed in tlie text. 



be injected into different input ports as shown in Fig. [5l 
Alternatively we could inject the two input fields C3,C4 
into the same port, but with opposite circular polariza- 
tions and replace the input beamsplitter with a polarizing 
beamsplitter. A half-wave plate is then inserted into the 
arm that does not interact with the atoms to rotate by 
TT the polarizations of 63 and 64 so that they are aligned 
with the polarizations of 0,3 and 0,4 respectively when 
they are recombined at a nonpolarizing output beam- 
splitter. With either version the mean difference current 
is zero because the optical fields are phase shifted with 
opposite signs and their contributions to the difference 
current cancel. 

We proceed by choosing A43 = A34 which gives 
K3 = K4. Solving the Hcisenberg equations for the 
light and atom evolution we find F°"* = R^(0at)F'", 

83°" = R5(^ph,3)S3'", S4™ = R^(^ph,4)S4"\ with the 
rotation angles 

e.t = 'i4[(iVph,4/2 + S'^l) - (iVph^3/2 + ^^3)] (9a) 
^Ph,3 = -~Ki{N,t/2~F'j') (9b) 
^ph.4 = -K4(^at/2-f Z'r). (9c) 

The rotation angles are a factor of 2 smaller than those 
of Eqs. ([2]) since now each frequency is assumed to only 
interact with one atomic ground state level. 

The initial conditions for the mean values of the pseu- 
dospin operators are (/'i") = N,t/2, (5^-^) = A^ph.s/S, 

{S™x) = ~A'ph,4/2, and all other components arc zero. 

The opposite signs of the values of (S™-^) and {S™x) 
precisely a result of the injection of the two fields through 
two different input ports to the interferometer. Note that 
in contrast to the single probe beam situation we now 
have (^at) = which implies that there is no differen- 
tial Stark shift of the clock transition due to the probe 
beams, provided they have equal mean photon numbers 
which wc will assume below. It follows immediately from 
Eq. ^ that the output variances of the light are 
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where now = |K4A'at-/Vph,4. The difTercnce photocur- 
rent at the output of the Mach-Zehnder is proportional to 



The expected value of the difference current is ) 
and the variance is 
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The variance includes a term proportional to 
K^A'at/-/Vpii,4 which is quadratic in the number of 
atoms. The reason the variance was strictly linear in 
the number of atoms for the single probe beam, but has 
a quadratic correction for two probe beams, can be seen 
by comparing Eq. ^Bp) with In the two probe 

beam case the rotation angles include an additional 
factor of iVat/2. Each probe is phase shifted proportional 
to the (fixed) number of atoms. However, the noise of 
the photocurrent difference depends on 6*32; — Si^x which 
has a coherent state variance. This gives a contribution 
to the measured photocurrent variance that is quadratic 
in the number of atoms. Thus the two- probe technique 
is suitable for observing atomic projection noise in the 
limit when A'at/A'ph.4 <C 1. Fortunately this limit is 
readily achievable in practice. 

From an experimental perspective the two input con- 
figuration has further advantages. Most importantly, the 
common mode noise e.g. arising from small displace- 
ments of the interferometer pathlength due to acoustics 
or vibrations, yields an opposite change in the output 
signal of the two probe fields. Hence, for equal power 
in the two input beams this noise will to first order be 
suppressed in the output photocurrent. In the config- 
uration where the probe fields enter the interferometer 
through two spatially separated input ports it is crucial 
that the fields have a very good spatial overlap in or- 
der that they interact with the atomic sample in exactly 
the same way. The configuration where the two fields 
enter the interferometer with orthogonal circular polar- 
izations on a beamsplitter facilitates this mode overlap 
more readily as the fields may be spatially overlapped in a 
polarization-maintaining fiber before the interferometer. 

In order to find the achievable spin squeezing we must 
again account for inelastic scattering. Before calculating 
the result for the Cs Dl line let us consider an idealized 
situation where the probe at ^^3 couples |3) |3') which 
only decays to |3) and the probe at uj^ couples |4) — > |4') 
which only decays to |4). In this situation all moments of 
the Fz operator are unchanged by photon scattering and 
the reduction in spin squeezing is only due to a reduction 



in the magnitude of the coherent spin state: {\Fy 
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77) I) . With the definition of spin squeezing relevant for 
Ramsey spectroscopy defined in Q we find 
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(1 - 77)2 1 + d-q 
where we have used — drj. From this equation one 
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FIG. 6: Optimized spin squeezing and scattering parameter 
as a function of optical depth d from Eq. (|14p . The inset 
shows the variation of with r) aX d = 100. With the same 
atomic parameters as in Fig. Uwe have again d = 22. Taking 
Afph = 9. X 10^ and A43/27r = 150 MHz we get 77 = 0.17 and 
2 - 0.34. 



may find the inelastic scattering rate that yields the high- 
est squeezing 



d-2 
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which for d ^ 1 gives r/o ^ 1/3. The corresponding 
maximal squeezing for large resonant optical densities is 
^2^i,j = 27/ (4c?). This 1/d scaling as opposed to the usual 
1/Vd is an attractive feature in the context of cold atomic 
samples with limited optical depth. 

Such an idealized situation is difficult to achieve in 
practice with available atomic level structures. One pos- 
sibility is to use the basis states \F,m = ±F) and the 
transitions \F, m = ±F) ^ \F' = F + l,m' = ±F ± 1). 
Such states may not be well suited for atomic clocks since 
the clock frequency is defined by Zceman shifts and will 
be both relatively small and linearly sensitive to magnetic 
field fluctuations. In the case of Cs we may alternatively 
use two-colour probing of the clock states |3), |4) with tt 
polarized light on the D2 line. With frequencies chosen 
such that / = 3 couples to /' = 2 and / = 4 couples to 
/' = 5 the effect of photon scattering is to populate states 
with ^ 0, but the value of / is not changed. Since 
the states with m/ ^ have a slightly different coupling 
strength to the light than the clock states some noise is 
added to the population difference measurement. It can 
be shown that for not too strong scattering such that 
we only need consider states with ruf = ±1 and optical 
depth not more than about 50 the relation ^^^jj^ = 27/ (Ad) 
is still a good estimate to an accuracy of about 20%. Al- 
though we cannot achieve the asymptotic 1/d scaling of 
the spin squeezing for very large d, we nevertheless re- 
tain the other advantages of two-colour probing discussed 
above. 

For a quantitative comparison with the one probe re- 
sults of Fig. [3] let us again consider two-colour probing of 
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FIG. 7; (color online) Arrangement for spin squeezing with 
amplitude modulated light. Solid lines are optical beam paths 
and dashed lines are electrical signals. 



the Cs Dl line as specified earlier in this section. Scatter- 
ing of a photon at ^3 couples |3) to |30), |31), |3 — 1), |41), 
or |4 — 1), with the kets labeled as |/m). Similarly scat- 
tering of an UJ4 photon couples |4) to |40), |41), |4 — 1), 
|31), or |3— 1). We denote the probabilities of these events 
by coefficients ij^^ for 103 and ■q'^^^^ for 014. For the Cs Dl 
line we find 7^30^ 



77/6, V3±i = ^/16, vi% = 577/48, 
and r][^ = 77/6, 77^^^ = 577/48, li^i = '?/16. The coef- 
ficients are normahzed so that J2f m vfm + '^/tn = 
in the discussion preceding Eq. ([7]) we assume that pop- 
ulation in the 777, 7^ states is removed from the system. 
Calculating as in Refs. [l^, [l^ we find 



,2 _ (1 - Iv)' 1 I 2 (1 - 

^ (1-77)2 l + d77 3'' (1-77)2 



(14) 



For d ^ 1 the optimum scattering probability scales as 
T]o ~ 1/%/^ and ^2^i,j ^ 1/Vd. Figure [6] shows the cal- 
culated spin squeezing as a function of optical depth. 
We see that for large d the spin variance is about twice 
smaller than for the single frequency probing of Fig. 31 



IV. SPIN SQUEEZING WITH AMPLITUDE 
MODULATED LIGHT 



An alternative approach to spin squeezing that does 
not rely on spatial interferometry is to measure the phase 
shifts between different frequency components of an am- 
plitude modulated light field. The basic scheme is shown 
in Fig. [71 By performing the phase sensitive measure- 
ments in the frequency domain there is no need for the 
atoms to be placed inside a mechanical resonator that 
has interferometric stability. This has the advantage that 
the relative phases of the different frequency components 
can be stabilized in a quiet environment that is physically 
separated from the atoms. 



We now have four optical frequencies interacting with 
the atoms where the two are symmetrically placed at 
a detuning of ±113/2 about uj^ and u)4± are symmetrically 
detuned by ±f24/2 about ^4. The central frequencies 0^3 
and UJ4 are chosen such that the magnitude of the phase 
shift is equal for light at the lower and upper sideband, 
typically implying that uj^ and LU4 arc set very near the 
resonance frequencies of the relevant atomic transitions. 
We emphasize that there is no light actually present at 
a'3,u;4 whereas the interaction strength and photon scat- 
tering rate can be adjusted as desired by variation of the 
detunings fl^jfli. 

The applicable continuous light operators are combi- 
nations of the sideband frequencies interacting with each 
transition. We introduce photon annihilation operators 
03+, 03-, 04+, a4_ and construct the operators 



Say 
S3Z 



-3 ("3+«3- - 0-3-^3+) tp 



3+^3+ - a3_a3_ 1 1 



•-pi 



and similarly with 3 replaced by 4. The photon num- 
ber operators once again are defined as -/Vph,3 = 

(j^3+o,3+ + o,3-0-3-^ tp, a-nd likewise for the F = A 
ground state. 

With the approximation of negligible photon scattering 
Niit,Fz, S3Z, and 84^ are constants of the motion. We 
can therefore integrate the Heisenberg equations for the 
light and atomic evolution to get F°"* = R^(^at)F'", 

83°" = R.z(^ph.3)S3'", 84°" = Rz(^ph,4)S4'", with the 
rotation angles 



Sat — K4S'^4 — ii-sS^^ 

9ph,3 = ;i3(^at/2 - Fi") 
9ph,4 = K4(7Vat/2 + Ff). 



(15a) 
(15b) 
(15c) 



where K3 and K4 are the common interaction strength 
of the respective pairs of sideband frequencies. For equal 
powers in the sidebands the initial conditions for the light 
operators are (S*^"^) = iVph,3/2, {S'^^) = iVph,4/2 with all 
other components equal to zero. As (S™^) = (S'^"^) = 
Eq. ([T5k ) displays that there is no change in the expec- 
tation values of the components of atomic spin F even 
when K3 7^ K4. In other words, the light shift is canceled 
by the combined influence of the two sidebands for each 
atomic level. 

However, for the atomic output variance we find 



where 



((AFr')' 



4 {klN, 



1 



ph,3 



ph,4. 



(16) 



(17) 
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The coupling constant is effectively four times bigger 
than for the Mach Zehnder scheme. The reason being 
that now all the light interacts with the atoms. 

For the light, the photo detector measures the com- 
bined power of all the involved fields. The detector does 
not respond to the high frequency interference between 
ijj-i and Lo^ so the photocurrcnt operator is proportional 
to 

i = (aj^ + aj_)(a3+ + a3_) + (aj_l_ + aj_)(a4+ + 04-) 

+TVph,4COs[fl4i + 4K4(7Vat/2 + F,)\. (18) 

The photocurrcnt is split in two and mixed with the 
local oscillators f23 and respectively to give two out- 
puts 

7Vph,3 sin[4K3(7Vat/2 - F^)] ~ ^ph,34K3(7Vat/2 - F,) = 03 
7Vph,4 sin[4K4(7Vat/2 + F^)\ ~ 7Vph,44K4(7Vat/2 + F^) = 6I4 

The difference of the measured phase angles is 

e - 4iVph,4^4(iVat/2 + F,)-4iVph,3^3(iVat/2-F,) 

= 2(K47Vph,4 - K3^'ph,3)^at + 4(K4^ph,4 + K3^ph,3)^.- 

(19) 

We can choose the coupling constants and intensities 
such that K4Afph,4 = K3iVph,3. With this choice 

Q = 2ki{N^^^i - ^PM7Vph,3)iVat + 8K4^ph.4^^z (20) 



ph,3 



which has a variance 



{AOf = K^iVat -I- 16(^:4)'iVph,4^at 



2iV, 



ph,4 



1 



2A^, 



ph,4 



(21) 



where is defined in Eq. pT]) . We must add to this the 
variance due to the shot noise of the light when no atoms 
are present and choosing for simplicity A^ph,3 = ^ph,4 
(implying that and f24 are chosen such that ^4 = ^3) 
the shot noise becomes iVph,3 -I- iVph,4 = 2A^ph,4, yielding 
an output variance 



{M f 2A^ph,4 



1 



^2 ( 1- 



at 



2A^ph, 



(22) 



Apart from a different numerical factor we find the same 
result as in Eq. (jlip for the two-probe interferometer. A 



quantum limited measurement of the phase angle 9 will 
project the atoms into a SSS with reduced variance of 
Fz and as in the two-probe case of Sec. IIIII there is an 
additional contribution to the variance which scales as 
-^at/^ph- The effects of photon scattering enter in the 
same way as in Sec. IIIII Thus this four-frequency ap- 
proach has the potential for good spin squeezing perfor- 
mance. The most challenging technical requirement is 
the need for quantum limited phase measurements at a 
frequency U.^ = ^4+ — 0^4- which must be at least a few 
times larger than the excited state radiative linewidth in 
order to keep the photon scattering rate sufficiently low. 
V. CONCLUSION 

We have analyzed the use of multiple probe frequen- 
cies for generation of spin squeezed atomic ensembles. 
Robust preparation of atomic spin squeezing requires 
suppression of technical noise and inhomogeneous light 
shifts, together with the freedom to choose the optical 
detuning in order to optimize the photon scattering rate 
for a given atomic sample size. QND interactions with 
a single probe frequency do not generally allow all of 
these requirements to be simultaneously met. We have 
shown here, using the Cs Dl line as a specific example, 
that the use of multiple frequencies with symmetrically 
chosen detunings can satisfy all of the above conditions 
simultaneously. Two possible configurations were pre- 
sented. The first uses a Mach-Zehndcr interferometer as 
in recent non-destructive measurements of the Cs clock 
transition [ol. [ll| . The second configuration uses only fre- 
quency domain instead of spatial interference which has 
the advantage of not requiring a mechanically stable in- 
terferometer. 

We have also pointed out that with multi-colour prob- 
ing on cycling transitions the spin squeezing variance 
scales as 1/c? as opposed to the 1 / \/7l scaling obtained 
with single colour probing. This type of multi-colour 
probing has been used in a recent demonstration of 
squeezing on the Cs clock transition [20j. 
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